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Abstract 

The Wheeler-DeWitt equation for the Bianchi Class A cosmological models is ex- 
pressed generally in terms of the second-order differential equation like the Klein-Gordon 
equation. To obtain the positive-definite probability density, a new method extending 
the Dirac-Square-Root formalism, which factorizes the Wheeler-DeWitt equation into 
the first-order differential equation using the Pauli matrices, is investigated. The solu- 
tions to the Dirac type equation thus obtained are expressed in terms of two-component 
spinor form. The probability density defined by the solution is positive-definite and 
there is a conserved current. The newly found spin-like degree of freedom causes the 
universe to go through an early quantum stage of evolution with agitated anisotropy- 
oscillation like Zitterbewegung. 

1 Introduction 

Arnowitt-Deser-Misner Q (ADM) reformulated the general relativity using the canonical for- 
malism, which makes a 3+1 split of the space-time metric, and introduced some constraints 
through the variational principle. According to the ADM approach, the Einstein equation 
is derived as the equation of dynamical evolution. The canonical formalism by the ADM 
approach is a system of dynamics including two constraints: one is the momentum con- 
straint and the other is the Hamiltonian constraint. The latter governs the time evolution 
of the system. According to Dirac g] the procedure of canonical quantization replaces the 
Hamiltonian constraint with the supplementary condition on the wave function which rep- 
resents the quantum state of space-time. The Hamiltonian constraint thus becomes the 
Wheeler-DeWitt equation, which is fundamental to the quantum gravity and determines 
the quantum state of space-time. The Wheeler-DeWitt equations for specific cosmological 
models like the Friedmann-Robertson- Walker model |3| or the Bianchi models [[![ |[ || have 
been investigated. 

The Wheeler-DeWitt equation for the Bianchi models is in general reduced to the second- 
order differential equation if the operator-ordering is neglected. The resultant Wheeler- 
DeWitt equation is similar to the Klein-Gordon type, which has a problem that the proba- 
bility density is not positive-definite. The Wheeler-DeWitt equation has also often this kind 
of problem that the probability density becomes negative. The negative values of probability 
density have been shown by Furusawa [Q through numerical calculation, who treated the 
Wheeler-DeWitt equation for the quantum mixmaster model using wave-packet solution. 
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To avoid the problem, the approach called Dirac-Square-Root formalism has been studied 
H El H' Hi ■ Since the Hamiltonian for the Bianchi type models is expressed in terms of the 
quadratic, the first-order differential equation can be derived by applying a procedure similar 
to the Dirac method. Mallctt ^ applied it to the Friedmann model coupled to the charged 
scalar field. Kim and Oh [Bl mainly discussed the Bianchi type-IX model. In their papers, 
the first-order equation included an unknown function dependent on the time-parameter is 
postulated and the wave function is required to be a solution of two-component column 
vector. Iterating the first-order equation thus postulated, the differential equation with 
respect to the unknown function is yielded so as to be consistent with the Wheeler-DeWitt 
equation for each component. The differential equation with respect to the time parameter is 
thus obtained. Indeed, the differential equation is reduced to Riccati or Bernoulli equation 
in case of the Friedmann model with charged scalar field or the Bianchi type-IX model, 
respectively. However, one sees that it is difficult to obtain the exact solutions in both cases. 

The purpose of this paper is to propose a new method, which factorizes the Wheeler- 
DeWitt equation into the first-order differential equation, using the Pauli matrices. The 
resultant differential equation thus obtained is similar to the Dirac equation. As a necessary 
condition for the Dirac factorization, it is postulated that the resultant Hamiltonian is "self- 
adjoint" . We apply this method to the the Wheeler-DeWitt equation of the Bianchi Class 
A vacuum models except for type-IX model. The solution of the Dirac type equation is 
written by two-component column vector. Using the Dirac type equation and the complex 
conjugate, we can define the probability density. Then, we get the equation of continuity and 
the conserved current. Here, we restrict ourselves to determine the form of Hamiltonian so as 
to make the probability density positive-definite. It is to be noted that our formalism cannot 
be applied to the Bianchi type-IX model, because the structure of anisotropic potential for 
the Bianchi type-IX model is different from that of the other Bianchi Class A models. 

The organization of this paper is as follows. In §2, we review the ADM canonical formu- 
lation for the vacuum Bianchi Class A models and introduce the Wheeler-DeWitt equation. 
In §3, we propose the new method in which the Wheeler-DeWitt equation is factorized into 
the first-order differential equation in terms of the Pauli matrices and apply the method to 
the Bianchi type-I, II, VIo, VIIo and VIII models. For each of them, the Dirac type equa- 
tion is derived and solved except for the Bianchi type- VIII model. Then, we assume a trail 
function which is expressed in terms of two-component spinor. We show that the probability 
density becomes positive-definite and there is a conserved current. As the time parameter 
approaches the limit of the singularity of the universe, the Wheeler-DeWitt equation for the 
Bianchi type-II, VIo, VIIo and VIII models are reduced to that of the Bianchi type-I model, 
so that the Dirac type equations must possess the solution for the Bianchi type-I model. In 
§4, we discuss the behavior for the Dirac type equation near the singularity in detail. It is 
shown that from the Heisenberg equation the quantum behavior of the anisotropic param- 
eters j3± make the quivering motion and the Hamiltonian system has a 'spin-like' degree 
of freedom. The 'spin-like' degree of freedom causes the universe to go necessarily through 
an early quantum stage of evolution with agitated anisotropy-oscillation. In §5, summary 
and discussion are made on the results of our formalism and comparisons are made with the 
Dirac-Square-Root formalism. Henceforward, we take that i, j and k runs from I to 3 and 
the unit c=%= I&tiG = I. 
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2 Canonical Quantization formulation for the Bianchi 
Class A models 

The Einstein-Hilbert action is given by |I(J 



1= J-"g w R*x, 



(4) (4) 

where g and R are the determinant of the space-time metric and the curvature, respec- 
tively; the indices (4) implies 4-geometry quantity. The metric for the Bianchi models is 
generally given by 

ds 2 = -N 2 dt 2 + e 2a e 2 ^ X l X J , (2) 

where N, a and are functions of t only: N is the lapse function; e 2a plays the role 
of a scale factor, so that the universe has the initial singularity as a — > — oo; /3y are an 
isotropic parameters of the universe, which gives a matrix form of traceless and x % are 1- 
forms. According to Misner 0], anisotropic parameters /3y are expressed in terms of two 
parameters j3+ and /3_ 

Pij = diag(/3+ + y/30_ , f3+ - Vl(3- , -2/3+). (3) 

The 1-forms obey 

dX* = \c} k x j A X \ (4) 

where Cj fc are the structure constants. 

The Bianchi models are characterized by the structure constants. The Bianchi models 
are classified into two groups from the structure constants, so-called Class A and B |TTJ. It 



should be noted that, as was first indicated by Hawking 12 , it cannot be reduced to the 
correct field equations by the ADM canonical formulation in the Bianchi Class B models 
except for the Bianchi type-V model [l4| . The reason is that the spatial divergence term 
does not vanish. Hence, we must discuss only the Bianchi Class A models because we use the 
canonical formulation. The characteristic of the Bianchi Class A models is that the structure 
constants obey the condition 

- 0. (5) 

We must therefore select the 1-forms so as to satisfy Eq. (|j). 

After the procedure for the ADM decomposition, we use the coordinates a, (3+, /?_ and 
the corresponding canonical momenta as p a ,p+,p-. From the variation of N, one obtains 

(3) 

the Hamiltonian constraint. The 3-geometry curvature R is obtained from the structure 
constants. According to Misner jjj, the relation of the curvature ( R to the anisotropic 
potential V(f3+, (3-) is introduced by 

(3) ii=|e 2 «(l-F(/3 + ,/3_)). (6) 

As a result, the Hamiltonian constraint becomes 

H=-p 2 a +p 2 + + P 2 _+e 4a (V-l). (7) 

On the canonical quantization procedure, we treat the operators p a , p+, p- in the 
Schrodinger representation: 

d - d . d 
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Because the Hamiltonian constraint is the first-class constraint, the quantized Hamiltonian 
constraint must be satisfied by the supplementary condition 



= 0, (9) 
f) 2 f) 2 f) 2 

H ^^-w + -m +e4a(v - 11 (10) 

where ^ = ^(a, /3+, /3_) is a state vector of the system, called by the wave function of the 
universe. This is the Wheeler-DeWitt equation for the Bianchi models. 

3 Dirac factorization for the Wheeler-DeWitt equation 

According to the Furusawa's discussion the probability density of Klein-Gordon type 

becomes negative for some values a. The negative values of the probability density are incon- 
sistent with the physical theories and it is the failure of the structure in the Wheeler-DeWitt 
equation. To get rid of the problem, we adopt the similar procedure that Dirac factorized 
the Klein-Gordon equation in order to obtain the positive-definite probability density. That 
is, we try to factorize the Wheeler-DeWitt equation into the first-order differential equation. 
We denote the first-order equation as Dirac type equation. As a necessary condition of Dirac 
factorization, it is postulated that the resultant Hamiltonian is "self- adjoint" , because one 
tries to think various ways of the factorization into the first-order equation and there is 
not any rules for the factorization. We suppose that the Wheeler-DeWitt equation for the 
Bianchi models is replaced by the equation of Dirac type 

H^> = (p a + axp+ + ct 2 ]3_ + a 3 e 2a v) = 0, (11) 



where v = y/V — 1 and the Pauli matrices <ii are expressed as 

1 \ / -i \ ( 1 



n < = ! i o )> ° 2= \ i o )> as v - 1 

We here rewrite the Dirac type equation (^lj) in the Schrodinger form 

d 

i-^> = [<riP+ + <y-it>- + <r 3 e 2a v] (12) 
The self-adjointness of the Hamiltonian is easily proved: 

7-0 = Pa + °iP+ + <?2P- + <y 3 e 2a v = H, (13) 

where the momentum operators are self-adjoint and commute with the Pauli matrices. 
With the help of relations (||) , the Dirac type equation (|Tl| ) is rewritten as 

<9* 5* d^> n 

~ i lT- ia i7^- ia2 J^ + a 3 e v ^ = °- ( 14 ) 
oa op + op- 

Taking the complex conjugate and multiplying vf^ to Eq. ( p^ ) from the left-hand side, we 
have 

i$t + i^ax— + i*V 2 — + &a 3 e 2a v$> = 0. (15) 
Oa op + op- 



Subtracting Eq. ( |l5| ) from Eq. ( |l4| ) multiplied ^ from the left-hand side, we obtain 

^ + _|_ ( ^ ) + _|_(^ )=0 . (16) 
Setting the probability density p and the current vectors j+ and j_ as 

p = j + = ^Vl*, j_ = * t (T 2 *, (17) 

we get the equation of continuity 

dp , #7+ 9j_ 



5a 9/3+ 0/3_ 



= 0. (18) 



Then, one sees that the probability density is positive-definite. We apply our formalism to 
solve the Bianchi type-I, II, VIo, VIIo and VIII models below. 

3.1 The Bianchi type-I model 

The metric of the Bianchi type-I model is given in Eq. (j^), in which the 1-forms are given 

by m 

X 1 = dx 1 , x 2 = dx 2 , X 3 = dx 3 . (19) 
The structure constants in the model obviously supply 

C l 3k = 0. (20) 

(3) 

It leads to the zero curvature of 3-dimcnsional hypersurfaces R\ = 0, which implies the 
spatial flat structure. From this, the anisotropic potential becomes V\ = 1. Then, the 
Wheeler-DeWitt equation is reduced to 

{-pl+f++V-)^ = Q. (21) 

We postulate here that the Dirac type equation for the Bianchi type-I model holds: 

Wi*i(o, P+,/3-) = (p a + aip+ + a 2 p-) * = 0. (22) 

The solutions to this equation have been already found JE], ^| . They are plane wave solutions 
given by two-component spinor form 



*i(a, = e i(k+0++k-0—^) ( </>i 

\ V2 



(23) 



where u> and k± are arbitrary real constants. Substituting the solution ( p^ ) into the Dirac 
type equation (p2|), we get the components 

, k + - ifc_ , , k + +ik- 

ipi = ip 2 and tp 2 = ipi • (24) 

to to 

It implies obviously I^Ail 2 = I V^l 2 and the relation among u and k± is 



u = ±Jkl + k 2 _. (25) 

Thus, the parameter u> corresponding to the frequency of the plane wave depends on the 
parameters k± corresponding to the wave number, which vary continuously. 



The probability density is 

ft (a,/3+,/3_) = = |^i| 2 + \H 2 = 2|^i| 2 , (26) 
which is positive-definite. The current vectors j± are respectively 

j + =^/>i+VrV>2 (27) 

and 

j- = iipstpi ~ iip*ip2, (28) 

which take some constants. Therefore, the equation of continuity dlq) holds and there is a 
conserved current. From the parameters (^) and the result of Eq. (|27D and (p8|), we interpret 
that the evolution of the quantized universe is that whether one of the three spatial axes is 
expanding and two axes are contracting or one axis is contracting and two axes are expanding. 
It should be noted that the quantity 

pi(a,j3+,p-)dj3 + dP- = 2|^| 2 / d(3+d[3„ (29) 



is not convergent and the wave function (|23|) is not the square integrable function. There- 
fore, it is impossible to normalize the wave function (|2^) without an appropriate boundary 
condition and we do not discuss the normalization constant in this paper. Here, the quantity 
(f29|) does not directly determine the probability, so that the probability must be regarded 



as a quantity proportional to this quantity (29). Moreover, since dp/da = 0, the probability 



density is independent of the time-parameter and quantity (|29|) is conserved. 

As mentioned above, it should be noted here that since the Wheeler-DeWitt equation 
and the Dirac type equation for the Bianchi type-II, VIo, VIIo and VIII models are reduced 
to that of the Bianchi type-I model at the singularity (a — > — oo), all models have the same 
property that the behavior near the singularity becomes that of the plane wave solution. We 
shall discuss the quantized Hamiltonian system of the Bianchi type-I model in §4. 

3.2 The Bianchi type-II model 

The metric of the Bianchi type-II model is expressed as Eq. (||) and the 1-forms are given 
by H 

X 1 = dx 1 - x 3 dx 2 , x 2 = dx 2 , x 3 = dx 3 . (30) 
Substituting the 1-forms ( |30|) into Eq. (|J), we get the structure constants 

C23 = — = 1. (31) 

Other components are zero. The curvature is 

(3) i? n = _I e 2«+4/3 + +473/3- 5 (32 ) 

and the anisotropic potential Vn (/3+ , (3- ) is obtained from the relation (^) as 

Vu-l = U iP++ ^ p -- (33) 
3 

Then, the Wheeler-DeWitt equation is expressed as 

-Pi +f++ f- + ±e 4a +^++ 4 ^- \y = 0. (34) 



G 



We postulate here that the factorized Hamiltonian is expressed as 



Hn =p a + a,p + + a 2 ^ + ?± e ^+W^- 



and that the Dirac type equation holds: 
With the help of the relation (0), it becomes 



■ggn 
da 



( J_ e 2Q+2/3 + +2\/3/3_ 

V3 



i d _ d 



\ 



d 



l_p2Q+2/3 + +2v / 3/3_ 



(35) 



(36) 



(37) 



Now we look for a solution to this equation. We take a trial function which gives the 
two-component spinor form and notices the form of the Dirac type equation as 



VPn = exp 



. 1 2a+2/3++2-v/3/3- 



6 



1p2 



(38) 



where the spinor components "01 and ip2 are arbitrary constants. Substituting the trail 
function (B8T) into Eq. (p7|), we get the relation of the spinor components 



t2 



-Ipl 



(39) 



Here we use the condition that the Dirac type equation ( |36[ ) for the Bianchi type-II model 
is reduced to that for the Bianchi type-I model as the limit of the singularity (a — > — oo). The 
trail function (38) would include such an asymptotic behavior. The trail function involving 
the plane wave part is expressed as 



%! = e i(k+P++k-P--va) exp 



2a+2/3 + +2 % /3 / 3- 



6 



-V>i 



(40) 



Although w, k± are arbitrary real constants, the coefficients to and k± are not free at all. 
Substituting the solution (Eo|) into Eq. we have 



— LU 

k + + ik- 



k + — ik- 

—U) 



Therefore, the coefficients must satisfy 

k, : 



-u), k- 



4>i 



0. 



0. 



(41) 



(42) 



Consequently, the solution is finally expressed in terms of the product of the amplitude and 
oscillation parts restricted by Eq. (f42|): 



*n = e - iu ^++^ exp 



1 2a+2|3 + +2v / 3,a- 
6 



-A 



The probability density is expressed, by using the wave function (43), as 



pa{a, p+, 0-) = tfjjtfn = 2|V-i| 2 exp 



1 2«+2(3 + +2i/3(3- 



(43) 



(44) 
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This is always positive-definite. One sees that this probability density goes to zero as f3 + + 
\/3/3- — > oo, while it approaches a certain constant as (3+ + \/3f3^ — ► — oo. The current 
vectors j± given by Eq. ( |l7| ) are 



3+ = -2|</>i| 2 exp 



_l_ 2a+2/3++2v / 3/3- 



and 



3- = 0. 



(45) 



(46) 



Hence, the equation of continuity (^8|) holds and there is the conserved current. Moreover, 
it is shown that the wave part propagates to the direction of (3 + axis only and it agrees with 
the result of Eq. 

3.3 The Bianchi type-VI model 

The metric of the Bianchi type-VIo model is expressed as Eq. (Q) and the 1-forms are given 
by p 



X 1 = cosh x 3 da; 1 — sinhx 3 da; 2 , \ 2 — ~ smn x 3 dx 1 + cosh a; 3 dx 2 , \ 3 = dx 3 . 
The structure constants are 

Q>3 = —C\ 2 = 1, Cgj = — C 2 3 = -1 
and the rest are zero. The anisotropic potential becomes 

4 



Vvio - 1 



(cosh(4v / 3/3_) 



1 



Then, the Wheeler-DeWitt equation for the Bianchi type-VIo model is expressed as 



2 + p\ + P- + ^e 4/3+ (cosh(4V3/?_) + l) 



* = 0. 



(47) 



(48) 



(49) 



(50) 



We postulate that the Dirac type equation for the Bianchi type-VIo model is expressed 



as 



^VI ^VIo 



Pa + axp+ + <T2P- + -^=a 3 e 2a+2 ^ (e 2 ^- + e - 2 ^ 



*vi = 0. (51) 



As in the case of the Bianchi type-II model, we assume a trial function of the form 



*vi = exp 



?/>! 
tp2 



(52) 



where the components ip\ and ip2 are arbitrary constants. Substituting the trial function 
(|52"|) into the Dirac type equation (E|), we get 



i>2 = -ipi 



(53) 



From the condition that the wave function ( |52|) should be reduced to the plane wave solution 
at the limit of the singularity, we obtain 



-oj, fc_ = 0. 



(54) 



The wave function involving the plane wave part is finally expressed as 



* VIo =e -**V>++") ex p 



_1 2c+2/3 + ( p 2V3f3- _ p -2V3P- 



The probability density becomes 
pvi (a,(3 + ,f3-) = H Io *vi Q = 2|^i| 2 exp 
which is positive-definite. The current vectors are 



1>l 



■2V3/3_ 



.7 + 



-2|^i| 2 exp 



2_ e 2 Q +2/3 + ^ e 2V3/3_ _ e -2\^/3- 



J- = o. 



(55) 



(56) 



(57) 



One sees that the equation of continuity holds and there is a conserved current. The plane 
wave part of the wave function (55) propagates to the direction of /L_ axis only. It is to be 
noted that since the wave function (p5|) and the probability density (|56| ) are both divergent 
as /?_ — > — oo, this solution is unphysical in such a region. 

3.4 The Bianchi type-VII model 

The metric of the Bianchi type-VIIo model is expressed as Eq. (^|) and the 1-forms are given 



by [15j 



X 1 = cosa^dx 2 + sinx 1 ^ 3 , % 2 = — sinx 1 dx 2 + cos x 1 dx 3 , y^=dx 1 . (58) 
The structure constants arc 



C23 



r 1 

32 



1, 

and the rest are zero. The anisotropic potential is 

4 



-a 



13 



-1 



VVlIo - 1 



3 



(cosh(4V3/3_ 



1 



Then, the Wheeler-DeWitt equation for the Bianchi type-VIIo model is expressed as 

4 



Pi + f+ + P- + 



fcosh(4\/3/?_) 



1 



* = 0. 



(59) 



(60) 



(61) 



Here, we postulate that the Dirac type equation with a factorized Hamiltonian is expressed 
as 



Hvio^VIIo = 



-Pa + &1P+ + &2P- 



V3 



CT 3 e 



2a+2/3+ ^2\/3/3_ _ e -2\/3/3_ 



*vn„=0. (62) 



As in the case of previous models, we can find out the wave function 

_l e 2a+2/3 + ^ e 2\/3/3_ + e -2\/3/3_ 



^ VIIo = e -MP++<x) exp 



The probability density becomes 

p VUo (a, /?+,/?_) = ^Ho^vnc = 2|^i| 2 exp \- 2 - e ^+ (e 2 ^- + 
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(63) 



(64) 



which is positive-definite. The current vectors have 

2 



j+ = 2|V>irexp 



-e 2 



3- = 0. (65) 



One sees that the equation of continuity holds and there is a conserved current. The plane 
wave part of the wave function (^3|) propagates to the direction of /3 + axis only. At the 
limit of /3_|_ — > —oo, the value of the probability density approaches a certain constant, while 
converges to zero as /3+ — > oo or /?_ — > ±oo. 

3.5 The Bianchi type- VIII model 

The metric of the Bianchi type- VIII model is expressed as Eq. (|^) and the 1-forms are given 
by HU 

X 1 — cosh a; 2 cos a; 3 da: 1 — sina; 3 c?x 2 , 

X 2 = cosh a; 2 sin a: 3 da; 1 + cos x 3 dx 2 , (66) 
X 3 = sinhx dx + dx 3 . 

The structure constants are 

C23 = — C32 = —1, C 31 = — C 2 3 = — 1, C 3 2 = — C 21 = 1 (67) 
and the rest are zero. The anisotropic potential becomes 

Will - 1 = ^ (e" 8/3+ - 2e 4/3 + + Ae- 2f3 + cosh(2%/3/3_) + 2e 4/3 + cosh(4V3/3_)) . (68) 



However, since the anisotropic potential (|68|) is a very complicated expression, we cannot 
present the Dirac type equation of the form as Eq. mA 



4 The behavior near the singularity 

It is interesting for us to discuss the behavior of the wave function near the singularity. 
Misner || discussed the Bianchi type- IX (Mixmaster) model which is represented by the 
classical dynamical system where an imaginary particle moves in the anisotropic potential 
of the approximated regular triangle and regarded the dynamical system as the billiard in 
the /?+-/?_ plane. In this way Misner's picture is very helpful to discuss the universe through 
the anisotropic parameters. 

The vacuum Bianchi type-I universe is so-called Kasner Universe. The evolution of 
the model in classical region presents an anisotropically expanding universe in which two 
axises are expanding and the other is contracting. While, in quantum region the quadratic 
Hamiltonian is derived by the ADM approach and the evolution of the quantized universe 
is presented by the Wheeler-DeWitt equation like Klein-Gordon type. There is the problem 
that the probability density can become negative by the Klein-Gordon approach. 

Let us consider to change the point of view to our formalism. Near the singularity 
(a — * —00), the Dirac type equations for the Bianchi Class A models are reduced to that of 
Bianchi type-I 

i-z- = (<riP+ + 02V- )^ = h'^f. (69) 
oa 

This dynamical system describes a massless imaginary particle moving in zero potential. 
The equations of motion in the Heisenberg representation for (3 + and p_ are 

= i h',$ + =io x p+,f3+ =o\ (70) 
da L J L J 
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and 



da 



h'J- 



1(72 



P-,0- 



02, 



(71) 



where the commutative relation \J3biPa] = ^ab (a, b = +, — ) is used. It is similar to the Zitter- 
bewegung |Iq| (a rapidly oscillating motion) of electron, so that the imaginary particle moves 
a trembling motions near the singularity of the universe. This behavior in quantum mechan- 
ics with factorized Hamiltonian is very different from that for the Wheeler-DeWitt equation. 
We emphasize here that this kind of behavior cannot be obtained from the Wheeler-DeWitt 
equation. 

Here, we set an orbital angular momentum rh' as 



rh' = (3+p- - (3-p+. 



(72) 



The Heisenberg equation of motion for the orbital angular momentum is 



We have, furthermore, 



dm' 
da 

da 3 
da 



h',rti 



°2P^ 



h' 



03 



= 2a 2 p+ - 2&iP- 



(73) 



(74) 



i<73 is a constant of motion. This result one can interpret by saying 

The Dirac type 



The quantity rh! 

a massless imaginary particle has a 'spin-like' angular momentum ^173 
equation (^) is a similar type to the Weyl equation, which has two-dimension of the spatial 
part here. We know that the Weyl equation describes the behavior of the massless neutrino 
with a half spin. On the analogy of this characteristic with the Weyl equation, the massless 
imaginary particle would have a half spin. Accordingly, we think that the Dirac type equation 
for the Bianchi type-I model describes the universe as a massless imaginary particle with 
'spin-like' degree of freedom. Thus, the 'spin-like' degrees of freedom cause the universe to go 
through an early quantum stage of evolution with agitated anisotropy-oscillation. However, 
the physical interpretation of the 'spin-like' degree of freedom is not comprehensible yet. In 
conclusion, we obtain that the Dirac type equation for the Bianchi type-I model shows the 
interesting results on quantum behavior, unlike Misner's negative conclusion that the status 
of the wave function remains classical near the initial singularity 0] . 



5 Conclusions and Discussion 

We have investigated the first-order equation using the Pauli matrices by factorizing the 
Wheeler-DeWitt equation for the vacuum Bianchi Class A models except for the Bianchi 
type-IX. We have derived the solution expressed in terms of two-component spinor on the 
Bianchi type-I, II, VIo and VIIo models. It has been shown that the probability density 
becomes positive-definite and the equation of continuity holds with a conserved current. 
However, we could not find the solution to the Bianchi type- VIII model because the form 
of anisotropic potential is very complicated expression. It is to be noted that we could not 
normalize the wave function and the probability density for each model. In the Bianchi type- 
II model, it is able to factorize the quadratic Hamiltonian of the Wheeler-DeWitt equation 
into self-adjoint Hamiltonians other than the Hamiltonian ([55]). However, we cannot extract 
physically interesting solutions to the Dirac type equations by such Hamiltonians. The 
similar results are obtained for the Bianchi type- VIo and type- VIIo models. 
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The Dirac type equations for the Bianchi type-II, VIo and VIIo model are reduced to that 
for the Bianchi type-I model near the singularity (a — > — oo). The quantized universe near 
the singularity is expressed by the plane wave solution with two-component spinor and the 
anisotropic parameters /3+ and /3_ make the motion of Zitterbewegung. The Zitterbewegung 
can be related to the 'spin-like' degree of freedom and thus the 'spin-like' degree of freedom 
causes the universe to go necessarily through an early quantum stage of evolution with 
agitated anisotropy-oscillation. The 'spin-like' degree of freedom which does not appear when 
we deal with the Wheeler-DeWitt equation for the Bianchi models may play an important 
role when we discuss the anisotropy of the universe in the initial stage. In such a sense, 
the Zitterbewegung is a significant subject how the big-bang or the inflation is related. To 
solve the problem of the initial singularity we must find out the physical interpretation of 
the 'spin-like' degree of freedom newly found by the Dirac factorization. 

As mentioned in §1, our method cannot be applied to the Bianchi type-IX model. Be- 
cause, the anisotropic potential of the Bianchi type-IX Vix is expressed in terms of 

Vjx = \e- 8t3 + - \e~ 2f3 + cosh 2^3/3- + \e^+ (cosh 4^3/3- - 1) + 1 > 0. (75) 

The quantity v 2 = Vix — 1 takes negative value as well as positive value here. It means that 
the quantity v — \fV\x. — 1 which appears in the factorized theory turns out to be imaginary 
as well as real. Therefore, we cannot factorize the quadratic Hamiltonian into a fixed form. 
Hence, we must not use the Pauli matrices but the Dirac matrices. To remedy this defect, 
we notice the fact that the anisotropic potential Vix itself is non-negative, and carry out the 
Dirac procedure of factorization with anticommuting 7^ (/x = 0, 1, 2, 3) and 75 matrices. 
One possibility is to make the factorization in the form 

(7V + 7 V + 1 2 P- + J 3 v'e 2a + 75 e 2Q ) * - 0, (76) 

where we set v' — %/^ix- However, as we have seen, because the anisotropic potential Vix is 
a very complicated expression, we could not have found out the solution to the Dirac type 



equation (76) with the factorized Hamiltonian. 

We compare our formalism with the Dirac-Square-Root (DSR) formalism || 9. On the 
DSR, it is postulated that the factorized Dirac type equation for the Bianchi type-I model 
is expressed as 

^-^w + -^w^ M+ ^- )=0 - (77) 

The Dirac type equation (|77]) then implies that the Wheeler-DeWitt equation for the Bianchi 
type-I model yields 

d 2 d 2 d 2 



da 2 dpi dpi 



tf(a,/3+,/3_) = 0, (78) 



for each component of the wave function. According to Mallett it is postulated that the 
Dirac type equation is expressed as 

i£ - Un± - U»± + iW(a)) *(«, MJ) = 0, (79) 

where VK(a) is an unknown function depended on a only. In Kim and Oh || paper, it has 
been applied the first-order equation ((79]) to the Bianchi type-IX model. 

However, it is not valid to apply the DSR formalism to the Wheeler-DeWitt equation 
including the anisotropic potential term like (F5I) and to factorize in the form of the Dirac 
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type equation (f79|). Because, it is inconsistent with the original Wheeler-DeWitt equation 
in the point that the differential operators with the anisotropic parameters /3+ and /3_ do 
not act on the unknown function W{a). The incorrect equation yields if one assume the 



equation ( 79 ) when the Wheeler-DeWitt equation has the anisotropic potential term included 
the anisotropic parameters. While, our formalism is valid to apply to the Bianchi Class A 
models because it is able to factorize the Wheeler-DeWitt equation included the anisotropic 
potential term into the Dirac type equation (|ll|). This is the extended DSR formalism. 

When we yield the original expression from the Dirac type equation discussed by our 
formalism and the iterating equation, we also obtain the extra terms which include the time 
parameter and the anisotropic parameters. Regarding the extra terms as a new constraint 
to the wave function and solving the constraint, we cannot get solutions which satisfy the 
Dirac type equation. As a consequence, we only discuss the Dirac type equation without 
imposing the constraint condition on the wave function. 
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